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The early steps of photosynthesis involve the excitation of reaction centres (RCs) and light-harvesting (LH) units by 
light. Historically, the electronic coherence across RCs and LH units has been consistently neglected as it does not play 
a significant role during the relatively slow energy transfer steps of the primary process. Here, we turn our attention 
to the initial absorption process, which occurs on much shorter timescales. We demonstrate that spatially extended but 
short-lived excitonic delocalization across RC-LH units plays a relevant role in general photosynthetic systems, as it 
causes a redistribution of direct absorption towards the charge separation unit. Using the complete core complex of 
Rhodospirillum rubrum, we verify experimentally an 80% increase in the direct optical absorption of the RC in situ as 
compared to isolated RCs. Numerical calculations reveal that similar enhancements can be expected for a wide variety 
of photosynthetic units in both plant and bacterial systems, suggesting that this mechanism is conserved across species 
and providing a clear new design principle for light-harvesting nanostructures 


Photosynthesis - the conversion of sunlight to chemical en¬ 
ergy - is fundamental for supporting life on our planet. De¬ 
spite its importance, the physical principles that underpin the 
primary steps of photosynthesis, from photon absorption to 
electronic charge separation, remain to be understood in full. 
Excitonic dynamics between tightly-packed pigments, such as 
within the RC or within the LH units, has been recognized to 
be of considerable importance for characterizing their individ¬ 
ual optical responses and for determinin^he associated time- 
scales for excitation energy transfer steps^®. For instance, the 
classical theory of resonance energy transfer between moder¬ 
ately coupled electronic excitations was originally formulated 
for transport amongst single pigmentJ^, and has subsequently 
required considerable modification for the treatment of exci¬ 
tonic delocalization among several pigments. The modified 
theor}® allows the estimation of energy transfer rates among 
RC pigments or between LH-bound pigments that belong to 
the LH unit which are in good agreement with the exper¬ 
imental observations. More recently, the study of coherent 
effects in these biologically relevant systems has increasingly 
attracted attention due to the observation of long-lasting os¬ 
cillatory signals measured with time-resolved techniqueP^ 
which have driven a wave of theoretical work able to success¬ 
fully explained these observation^®^. These studies have 
neverthelss confirmed that, even when assuming the electronic 
coherence between RC and LH units to survive hundreds of 
femtoseconds, its effect on energy transfer processes between 
LH and RC, which take place on time scales of tens of pi¬ 
coseconds, is rnino]!^. As a consequence, quantum mechan¬ 
ical excitonic delocalization across RC and LH units has so 
far been consistently ignored. However, we will show that 
for processes that are taking place on much shorter timescales 
than energy transfer, excitonic delocalisation among RC and 


LH pigments may become relevant. 

The principles of directed absorption The fundamental 
principle that motivates both theory and experiment in this 
work is illustrated in Fig. lA. Light absorption in photosyn¬ 
thetic complexes occurs predominantly at LH antenna com¬ 
plexes, composed of many pigments, followed by excitation 
energy transfer to the RC, containing far fewer pigments. The 
process of photon absorption is described by the absorption 
spectrum, which quantifies the probability of exciting specific 
transitions and ultimately, contains signatures of the excitons 
in an absorbing unit at early times. The excitonic delocal¬ 
ization across RC and LH pigments exists for times deter¬ 
mined by the decoherence processes (~ 100 fs) which are 
much longer than the few optical cycles (~ 10 fs) required to 
imprint its effects on the absorption spectrum. Consequently, 
even for the case in which excitonic energy transfer (EFT) is 
well-described as a classical transport process, quantum me¬ 
chanical delocalization of excitons across the entire complex, 
antenna and RC (see l.h.s. of Fig. 1 A), must be accounted for 
in an accurate description of the absorption process. After a 
few hundred femtoseconds, environment-induced dephasing 
will degrade exctionic coherence (see r.h.s. of Fig. lA) and 
the subsequent EFT can be accurately described by incoherent 
classical hoppinj^. We will show that even if EFT is incoher¬ 
ent, the initial quantum mechanical delocalization in the ab¬ 
sorption process can lead, remarkably, to a significant redistri¬ 
bution of absorption intensity from antenna to the RC, thereby 
increasing the proportion of excitation which produces charge 
separation without the requirement of further transfer steps. 

Our model system for exemplifying this mechanism theo¬ 
retically and verifying it experimentally is the core complex of 
the purple bacterium R. rubrun^. This organism is a fresh¬ 
water bacterium, which is found in shallow, turbid pondss un- 
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der conditions of low light intensity and low oxygen tension. 
Under these conditions the overall light intensity reaching the 
purple bacterial population at about 2 m depth is extremely 
low, as indicated by recent in situ spectral measurement^. 
We estimate that the the incident light intensities could be as 
low as 60 photons.nm^^.ms which may be severely lim¬ 
iting, considering that the time for one turnover of photoin- 
duced cyclic electron transport is about 1 mJ^. Thus, R. 
rubrum is an interesting organism for discovering design prin¬ 
ciples which have arisen as a response to the strong evolution¬ 
ary pressure occurring under these low light conditions. 

The core complex possesses a photosynthetic RC, con¬ 
taining a “special pair” (P) of bacteriochlorophyll (BChl) 
molecules which mediate the primary process of light-induced 
charge separation, surrounded by an LHl complex contain¬ 
ing 32 BChls arranged in a protein with ring-like geometry 
around the RC (Fig. lA). The absoprtion spectrum A(/i) of 
the P pigments peaks around 870 nm, whereas the LHl ex¬ 
hibits a single absorption band at 880 nm, denoted the P870 
band and B880 band, respectively (see Fig. IB and Supple¬ 
mentary Information for spectra). The single LHl band arises 
due to the circular arrangement of the 32 BChls with the Qy 
transition dipoles lying almost tangentially to the LHl cir¬ 
cumference. The spectral bands of the P pigments or LHl 
correspond to transitions from the electronic ground state to 
the excited states |P870) and |B880), and their strengths are 
determined by the effective dipole moments Dpgjo and Desso, 
respectively, according to A(/i) oc (see Supplementary 
Information). 


RESULTS 

One tantalizing and experimentally verihable consequence 
of our full Monte-Carlo numerical simulations is a signihcant 
increase of the direct absorption of the P pigments when in the 
core complex (Fig. IB). The full numerical simulations take 
into account all the pigments which compose the core com¬ 
plex, subject to variations in energies and coupling strengths 
induced by the highly dynamical protein scaffold (details in 
the Supplementary Information). However, the key result can 
be understood by a simple model with two effective excited 
states |P870) and |B880) coherently coupled with strength 
V through the Coulomb interaction between their transition 
dipoles, as described by a Hamiltonian of the form 

■K = A£ |P870> <P870| -h V (|P870) <B880| -h |B880) <P870|) , 

( 1 ) 

where AE - Ep^jo - Fbsso > 0 is the inter-band energy differ¬ 
ence between the two bands. The eigenstates la) of the system 
Hamiltonian in Eq. (1) are delocalised across both the RC and 
the LHl and lead to dipole moments Z5', which become linear 
combinations of Dpg^o and Db&so- To lowest order in (VIAE), 
valid when V « AE, (full expressions in Supplementary In¬ 
formation) this redistribution of dipole strength can be written 

|T>p87oI^ “ |77p87ol^ + 2 Dpsjo ■ Dsm (^) + — (2) 
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Figure 1. Evolution of coherence and the relevance of differen¬ 
tiated time scales. (A) Schematic representation of the R. rubrum 
RC-LHl bacteriochlorophyll (BChl) pigments. The positions of the 
LHl and RC BChls have been obtained from a molecular model (see 
Supplementary Information). After a few hundred femtoseconds, 
environment-induced dephasing will degrade the quantum mechani¬ 
cal delocalization between the antenna and the RC (red shading), and 
transforms into a statistical mixture of populations residing either on 
the antenna or the RC (blue shading), meaning that absorption is co¬ 
herent while transfer is not. (B) Experimental and theoretical absorp¬ 
tion spectrum of the R. rubrum RC. Extended delocalization across 
the RC-LHl leads to an 80 % increase in absorption intensity of the 
P870 band, as illustrated by the area (yellow) between the spectra for 
the core complex and for purified RC. 


The coherent nature of the coupling is manifested in two mea¬ 
surable optical signatures: spectral shifts and absorption in¬ 
tensity redistributions, which scale as (VjAEY and VIAE, re¬ 
spectively, to lowest order. Hence, in the core complex, where 
coupling among RC and LHl satisfies V <Si AE, excitonic 
coupling may lead to a measurable absorption intensity re¬ 
distribution while leaving spectral positions almost unaffected 
due to the quadratic dependence on (VIAE). 

According to equation (2), absorption of the P870 band is 
enhanced when the ratio (Dp%iQ ■ i3B&s.o)(VIAE) is positive. 
Under these conditions, the absorption of the B880 band de¬ 
creases in equal measure such that the total absorption inten¬ 
sity remains constant. The geometrical details of the core 
complex imply that (Dpsio-Db&so) F > 0, a direct consequence 
of the almost tangential arrangement of LHl BChl transition 
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Figure 2. Comparison of the absorption redistribution between 
experiment and theory. (A) Experimentally determined difference 
spectra of pigment-protein complexes under different oxidation con¬ 
ditions. (B) Theoretical prediction of redox difference spectra shown 
in (A). Yellow shading indicates a redistribution of absorption which 
results from increased P870 absorption and a concomitant decrease 
in that of the B880. 

dipoled^2122l Therefore, Afs > 0 provides an energy landscape 
that inhibits energy transfer funneling to the RC, enhancing 
P870 band absorption at the expense of slowing the LHl—> P 
transfer step. 

Experimental confirmation of directed absorption The 
experimental confirmation of the predicted enhancement of 
the P870 absorption due to quantum excitonic delocalization 
is complicated by the fact that the P870 absorption spectrum 
overlaps with the much stronger optical response of the LHl 
B880 band. We carried out two independent procedures for 
the oxidation reaction, P —> P^, namely, optical excitation or 
stoichiometric addition- in the dark- of the oxidising agent 
potassium ferricyanide as explained in detail in the Supple¬ 
mentary Information. The oxidation reaction moves the P870 
resonance to 1260 nm, far into the infrared region, and there¬ 
fore effectively suppresses the resonant coherent interaction 
between P and LHl pigment^. Subtracting core complex 
spectra of oxidised P^ from those of neutral P filters out the 
overwhelming B880 component and highlights the changes 
in the absorbing bands due to the coherent interaction V. 
With both procedures we observed difference spectra from 


core complexes and purified RCs (Fig. 2A) which are in good 
agreement with theoretical calculations (Fig. 2B), and consis¬ 
tent with an 80 % increase in P870 absorption in the former 
over the latter (Fig. IB). 

The enhanced absorption of the P870 band leads to the 
excitation of a state that is partially delocalised over P and 
LHl pigments, so it is not immediately clear that increased 
absorption at 870 nm implies the increased excitation of the 
P-pigments, which is desirable as it places excitations di¬ 
rectly where charge separation takes place. Even though the 
BChls have several optical transitions besides the Qy transi¬ 
tion, under physiological conditions, the Qy transition(s) is the 
most relevant in R. rubrum. Indeed, recently published com¬ 
plete spectra from different water depths in oceans and lakeJ^ 
showed that in coastal waters and peat lakes blue light is ab¬ 
sent, and that the spectra of these ecosystems are dominated 
by red and near-IR light. As such, under incoherent broadband 
illumination covering the entire near-IR spectral range, the 
prominent excitation concerns the Qy transitions. Thus, un¬ 
der physiological conditions, the excitation population of the 
special pair, Pp(Q), is determined by the absorbtion from the 
Qy transitions, accounted for by a statep(O) oc 2^, \D'J^ \a) {a\. 
The requirement that many harvesting pigments should sur¬ 
round fewer charge separation pigments in order to increase 
the excitation of charge separation pigments, is illustrated by 
an analytical expression for Pp{Q) - Tr {p(0) IF) (FI) in the 
absence of protein-induced inhomogeneities 

■) 

(3) 

where M represents the number of excitonically-coupled har¬ 
vesting units encircling the N pigments that initiate charge 
separation in the RC (see Supplementary Information). In 
the particular case of the core complex these correspond to 
M - 32 pigments for the LHl and N - 2 pigments for 
the special pair. Full numerical simulations of the complete 
core complex result in a two-fold average increase of popu¬ 
lation PpiO) due to extended delocalization over antenna pig¬ 
ments, and it provides the initial condition for dynamical sim¬ 
ulations of EET (see Eig. 3). The probability of charge sep¬ 
aration, Pp+{t), increases as extended delocalization focuses 
the excitation into the special pair, allowing the excitation to 
avoid the slow transfer step from LHl to RC. This increase in 
charge separation provides a quantitative link between a co¬ 
herent quantum effect and enhanced functionality. 

It is important to note that the physiological environment 
of light-harvesting complexes is highly disordered and dy¬ 
namic, leading to variations in the excitonic energy splitting 
AE. Upon averaging the complete expressions that lead to 
perturbative equations (2) and (3) over appropriate distribu¬ 
tions (see Supplementary Information), Eigure 4 shows that 
excitation population is redistributed from antenna to F pig¬ 
ments for all realistic levels of disorder. Moreover, a maxi¬ 
mum occurs for intermediate values of disorder, equal to half 
the excitonic energy splitting, AE 12. Remarkably, the param¬ 
eters which are consistent with pump-probe studies in RC^ 
and reproduce our experimental dark-minus-light spectra (see 
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Figure 3. Absorption redistribution enhances charge separation. (A) Transport dynamics in R. rubrum can be accurately modeled using 
classical rate equations (see Supplementary Information). (B) Both the initial excitation of F-pigments (Pp(t) (dashed lines)) and the resulting 
process of charge separation, Ppt (t) (solid lines), are increased when absorption is delocalized across the core complex (red) as compared to 
absorption occuring either on the RC or LHl (blue). The enhancement is relevant on a timescale comparable to the charge separation process 
(3 ps, grey line), well before EET from LHl to RC takes place (~40 ps). The initial excitation of R-pigments increases from 2/34 to 4/34 (an 
increase of 100%), and results in a ~ 1% increase in photosynthetic efficiency at long times. 


Supplementary Information) optimize the population redistri¬ 
bution to the P pigments. Accordingly, neglecting disorder 
underestimates the 80% enhancement in the RSTO direct op¬ 
tical absorption, presented in Figure 1, due to this non-trivial 
enhancement arising from protein-induced inhomogeneities. 

Generality of directed absorption in nature The redirec¬ 
tion of excitation towards the special pair due to extended 
quantum delocalization can be generalized to other organisms 
by theoretical analysis. When environmentally-induced dis¬ 
order is included, numerical simulations predict similar re¬ 
distributions of excitation population in other purple bacteria 
species {N - 2, M - 2% - 32) which show structural devi¬ 
ations from a closed circular ring. In photosystem 1 (PSl) 
(see Fig. 4) of higher plants (N - 2, M - 83), our analy¬ 
sis indicates a 14 % increase of the excitation population of 
the P-pigments due to excitonic delocalization over antenna 
pigments. For the photosystem 2 (PS2) we predict a 12 % in¬ 
crease in P-excitation for both the monomer {N -2, M - 35) 
and natural dimeric structure {N - A, M - 70) (see Supple¬ 
mentary Information for further details). 


DISCUSSION 

Natural optimization of R. rubrum core complex It is in¬ 
teresting to examine the extent to which the design of the 
LHl-RC unit is optimized with respect to the population re¬ 
distribution due to excitonic delocalization. Figure 4 analyses 
the situation of ring-like LH units and shows that the opti¬ 
mal population redistribution occurs when the RC contains a 


small number of pigments, N, able to initiate charge separa¬ 
tion, while the LH unit contains a few tens of pigments, M, 
i.e. /V «c /W, as expected from Eq. (3). A greater number M 
of harvesting pigments naturally increases the distance of the 
LH to the RC pigments, reducing the coherent coupling, V, 
and therefore weakening the population redistribution, also il¬ 
lustrated by Eq. (3). Photosynthetic complexes typically con¬ 
tain N - 2 pigments which promote charge separation. In 
this case, the maximal enhancement for direct excitation is 
achieved with M ~ 20 pigments in the LH unit, closely re¬ 
sembling the LHl structure. 

Design principles: foccusing absorption in artificial de¬ 
vices A quantitative understanding of light harvesting in the 
LHl-RC complex illuminates clear design principles for de¬ 
vices engineered to direct light collection. These principles 
could prevent photobleaching and lengthen the lifetime of har¬ 
vesting devices by directing absorption to stable target pig¬ 
ments, in addition to increasing the efficiency from photon 
absorption to charge separation. By tuning the resonance con¬ 
dition AE in a ring-like geometry, a 400% increase in the ab¬ 
sorption of a few target pigments can be achieved. Moreover, 
in devices where dissipation is fast (~ 15 ps as opposed to 250 
ps in LHlEl), these configurations will lead to a 54% gain in 
achieving charge separation following photon absorption (see 
Supplementary Information). 

In conclusion, theoretical considerations have led us to 
identify a new design principle which exploits quantum de- 
localization to redirect light absorption to the RC in several 
photosynthetic systems. Given the generality of the design 
principle described here, it is tempting to equate this princi- 
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Figure 4. Inhomogeneous broadening and clustering of LH and P pigments as design principles for optimal performance. Exci¬ 
tation population is redistributed from antenna to P pigments (A) for levels of disorder up to the full spectral linewidth (sld ~ 200-300 
cm“*). The redistribution is, remarkably, noise-assisted, as the optimal redistribution occurs in R. rubrum for intermediate levels of disorder 

( yjstdp + std\u^ = A£/2 = 65 cm“'. PS 1 (B) and PS2 (see Supplementary Information) exhibit clustering of many LH pigments around fewer 
RC pigments (bound in dark blue). The excitation population of the P-pigments (dark blue arrows) increases when extended delocalization 
over all antenna pigments (bound in green) is taken into account. Error bound corresponds to the width of distribution when averaging over the 
effects of environmental disorder. Using a simplified model for a general circular aggregate (C), the percent increase of the initial excitation 
of the P-pigments is calculated (D) as a function of the number N of P-pigments and number M of LH pigments, and (E) as a function of the 
number M of LH pigments fixing N = 2. The enhancement of direct P-pigment excitation requires N <K M and is maximized for M ~ 20. 


pie with increased evolutionary fitness of certain biological 
systemJ^. We have verified this principle experimentally in 
the core complexes of R. rubrum, and our detailed theoretical 
analysis illuminates clear guidelines for the design of light¬ 
harvesting devices that direct photon absorption. This new 
quantum design principle provides a quantitative link between 
structure and function, broadening our understanding of the 
role of quantum mechanics in light harvesting. 
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ditional support was provided by the National Science Eoun- 
dation through PPC@JQI and the National Science Eounda- 
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Appendix A: Materials and Methods 


1. LH1 structure and Hamiltonian 


In this study we have employed a model of the R. rubrum 
LHl where 32 BChl molecules, bound to 16 a and /3 polypep¬ 
tides as QfjS(BChl )2 subunits, are arranged in a ring-like ge¬ 
ometry with Ci6 symmetry, see Fig. SI. The model 
was constructed in the group of one of the present authors 
(R. GhosP) in a collaboration with the group of Klaus 
Schulten (Beckmann Institute, University of Illinois, Urbana- 
Champaign) in 2002, and was based upon both extensive 
biochemical and biophysical data, low-resolution projection 
structures obta ined using cryoelectron microscopy (cryoEM) 
of 2D crystalJ^I^, and homology modelling to the known 
X-ray structures of the LH2 complexes from Rhodospiril- 
lum molischianun^ and Rhodopseudomomas acidophi 7^3 as 
well as a further homology model of the LHl structure from 
Rhodobacter sphaeroideM The data from this model, as well 
as that for the complete R. rubrum photosynthetic unit (PSU) 
(containing the RC) has been used subsequently for biophysi¬ 
cal studies of this complexP® 


The Ci6 symmetry naturally arises from the energy mini¬ 
mization procedure needed to construct the model, and is one 
of the controversial aspects since it contradicts AFM and cry¬ 
oEM data obtained for LHl complexes of other organisms^EI 
However, fluorescence polarization data obtained from both 
purihed R. rubrum LHl and LHl-RC complexes, respectively, 
reconstituted and diluted into bilayer membraneJ^, have con¬ 
sistently yielded results consistent with Cie symmetry, so we 
feel justifed to employ this assumption for the calculations 
made here. Very recently, the X-ray structure of a ther¬ 
mostable LHl-RC complex from Thermochromatium tepidum 
at near-atomic resolution of 3.0 A has been publishecp2| This 
LHl-RC complex is also comprised of 16 Q'yS(BChl )2 units 
which form a closed ring, und thus strongly resembles the R. 
rubrum PSU model employed here. However, the T. tepidum 
complex also shows overvable ellipticity. Nevertheless, the 
distances and geometries of the BChl pigments do not differ 
strongly from our R. rubrum model. Furthermore, the crys¬ 
tals were grown in detergent, which may have a destabilizing 
effect upon the geometry of the LHl-RC complex. 


In this study, we have employed the assumption of Cie sym¬ 
metry for the structure of the LHl. Its pigments are identihed 
by indices (m, 0, where m 6 {-M/4 + 1,..., M/4) labels the 
(dimerized) unit cell (M = 32), and ( e {1,2) labels the pig¬ 
ment within the a/3 heterodimer. Low energy spectral features 
can be characterized using the Qy optically-induced transition 
dipole moments of the BChl, which lie along a line joining 
the nitrogens of unreduced pyrrole rings I and IlPi! A gen¬ 
eral parameterization of the LHl BChl positions, and the 


A 




.1 


I 


. X 

I Acc BChl P BPh 1 
A , BPh P Acc BChl * 

X f 


B 


BChl' Zf Acc BChr*\^ 
BPh|, ' *BPh^^^ 


Figure 5. (A), (B) The position of the BChl and bacteriopheophytin 
pigments in the atomic model of the R. rubrum LHl-RC complex 
used here. The orientations of the Qy transition moments of the pig¬ 
ments are shown with arrows. The special pair (P), accessory BChls 
(Bl and Bm) and bacteriopheophytins {Hl and Hm) of the RC are also 
indicated. The a and (3 protein subunits are rendered translucent. 


corresponding Qy transition dipole moments, dm,^, is: 


dm,( — d 


^ mX 


-(-l)^sin(^m-H ^(-l)^-HA 7 f)cos< 
(-l)fcos(^OT-H |(-l)f-I-A7f)cos^ 
sin (p^ 

rfCos(^m-H ^(-l)f)' 
r^sin(|m+ |(-l)f) 

{-lYz/l 


(Al) 


An explanation of structural parameters and the values used in 
the full numerical study are summarized in Table S1. A corre¬ 
spondence between microscopic structure and optical proper¬ 
ties is possible by a theoretical analysis using a Hamiltonian 
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coupling Qy transition dipoles. We consider the Hamiltonian: 

m,^ 

(A2) 

where e„,,f is the site energy of BChl (m, and ^ is the in¬ 
teraction energy between pigments {m,() and {m 'When 
the details of the charge distribution of excited states are of 
no importance, j"^’^ can be calculated using the point-dipole 
approximation: 


necessary parameters are given in the Supplementary Table 
S 6 . The Hamiltonian of equation |A2| is best examined in the 
Fourier basis 


\K0 


1 

^[MJ2 


\m,0 


(A 6 ) 


where k e {-M/4 + 1,..., M/4). In this basis, the states de¬ 
couple for different values of k 

{k,^\H\k',n^k.('(k)Sk,k' (A7) 

Lu'(k) = 2 ^ I'”’ 


J 


-m'_ 




AnK 


in, ■ L,' 3 [L, ■ AC,/') [L',' ■ AC'/') 


Ar'/I 

m,4 


\ 


where hr"'f - 


(A3) 

rm, - and k is the relative permittiv¬ 
ity. The sub-nanometer distance between neighbouring chro- 
mophores implies that the nearest-neighbour couplings de¬ 
pend on the geometry of the electronic wave-function of each 
chromophore. These couplings can be fitted using fluores¬ 
cence anisotropy measurements, resulting in Qi = 600 cm“^ 
and Q 2 - 377 cm ' for the intra- and inter-dimer couplings, 
respective!)®. 


Diagonalization of equation A2 leads to excitons |a) = 
Yjm, c“ ^ \m, ^), where a is a collection of quantum numbers 
which describe the state. The optical response depends on the 
dipole moment of the exciton states: 




(A4) 


and is calculated by means of the Fermi’s Golden Rule-type 
expression 




(A5) 


where homogeneous broadening is taken into account with 
a lineshape function, /. The brackets correspond to ensem¬ 
ble averages that take into account inhomogenoeus broaden¬ 
ing over distributions of fluctuations in the site-energies and 
nearest-neighbour couplings. 


2. Analytical Results 
a. Dipole Redistribution 

The exciton states of circular aggregates are well-knowrP^, 
and modeling the LHl complex with a Cie symmetry allows 
simple analytical expressions for noiseless spectral observ¬ 
ables. To make these expressions transparent, we consider a 
simplified model whereby the transition dipoles sit on a two- 
dimensional ring of radius R, the BChl are degenerate in en¬ 
ergy e and share a common dipole moment strength cfi, and 
only nearest-neighbour couplings are taken into account. The 


The diagonalisation of the 2x2 Hamiltonians restricted to 
each k-subspace, L(k), is accomplished with the unitary: 


cos0(k) e®W/2gin0(k) \ 

U(k) - j (Ay) 


where ^(k) = arg(Loi(^)) and ®{k) = ^ arctan re¬ 

flect the “amount of dimerisation”: fl)(k) captures the differ¬ 
ence between intra- and inter-dimer coupling, and 0 (k) cap¬ 
tures the site energy differences. 

Considering degenerate dimers and taking into account 
only nearest neighbour coupling, the matrix L takes on a par¬ 
ticularly simple form: 




e 01 + Qie 

01 + 02^“®^ e 




(AlO) 


where e, 0i, 02 are site energies, and intra-dimer and inter¬ 
dimer couplings, respectively. Accordingly: 


m) = ^ 


O/k) = arctan 


02 sin ^k 
01 H- 02 cos ^k 


(All) 


The Hamiltonian is diagonalized by the exciton wavefunc- 
tions |k, v), leading to dipole moments Z5^y and energies 
E(k, v): 


\k,v) = 


1 

^[MJ2 


y]e's'='"t/(k)f,Jm,/'> 

m,^ 


Dk,i - Sk,±id COS 

Dk,2 - dk,±id Vm /2 sin I 


<F(1) -I- y\ e 


V2 


^(i) + e 


V2 


V2 

v? 


An, 


£(k,v) = e + (-l)''| 0 i+ 02 Cos—k|, 


(A12) 

(A13) 

(A14) 

(A15) 


where y is the intra-dimer angle and Ay is the tilt of the dipoles 
away from the tangent, as described in Table S1. The circular 
symmetry is manifest in the two-fold degeneracy Fj. ^ = E_k,v 
and the dipole selection rule D^ y cx 6k,±\, and the dimerisation 
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in the splitting by 2(Qi - Q 2 ) of the energies into two bands, 
denoted by v. 

The concentration of dipole strength in the upper and lower 
bands is a function of il)(l) + y, which reflects the amount 
of dimerisation. Increasing dimerisation, Q\ > Q 2 , leads to 
0 ( 1 ) + 7 « 1 and, consequently, more dipole strength in the 
lower band. For the LHl, more than 97% of the total absorb¬ 
ing strength is in the low-lying \k - ± 1 , 1 ), which is respon¬ 
sible for the strong response of the LHl at £’(+1,1) = 880 
nm. 

The functional advantage of this dipole strength concentra¬ 
tion becomes clear after calculating the coupling between the 
LHl and P-pigments. The P870 band corresponds to absorp¬ 
tion by the state |P870) = (|P1) - |P2)) since the dipole 

moments of the P-pigments are nearly anti-parallel. This 
bright special pair state is modeled by; 


where N — 2 and /? parameterizes the dipole moment of |P1), 
approximately lying in the plane of the ring. The coupling 
between LHl excitons and P-pigments, is then, after some 
algebra using equation | A3 1 


VmN cP- 

y,,i =<P870|7|^,l) = 5,,±i^^^cos 
X (1 + sin Ay) 


0(1)-Hy 


(A17) 


y,,2 = <P870|7|^,2) = 5,,±i- 


mN 


cP . / 0 (l)-t-y 
^ sm - 


2 P3 ^ 

X (+/)e±'03-Ar) ^ (a18) 


where R is the radius of the LHl ring. These expressions il¬ 
lustrate a selection rule, Vk.v 6k,±\, for the coupling between 
LHl excitons and the P-pigments, and show that increased 
dimerisation enhances the coupling between P870 and the 
near-resonant |+1,1), boosting Forster rates. Since |P) only 
effectively couples to \k - +\, 1 ), we need only three-levels in 
order to characterize RC-LHl interactions. In what follows 
we drop the index v and consider the exciton states and their 
associated dipole moments; 


|A: = +1) = - \ - g '®(i)/2 \ffi,2}) 

Vm V ^ ’ 

(A19) 
(A20) 


4 =±i = d 4Mn cos j 


/— i \ 

,V2 


V a/2 / 


The symmetry of the problem leads to a selection rule cou¬ 
pling only |P870) and \k - +1), so the three-level Hamiltonian 
required to capture RC-LHl interactions is; 


H = 


A£ Fi y_A 


y* 


0 

0 


(A21) 


A17 


This 


where AP s 130 cm ' and Vi = Fi j of equation 
three-level system can be reduced to an effective two-level 


system by choosing a suitable basis for the degenerate k - ±l 
subspace. This basis is found to be 

|P880> = (e”''* = 1) H- e''* = -1)) , (A22) 

V 2 ^ ’ 

|0> = — = 1) - P'’' \k = -l}) (A23) 

V2 ^ ’ 

where i/r = arg(yi). Note that |£)b 88 oI^ = y- In that basis the 
Hamiltonian becomes 


H = 


'AE V 0) 

y 0 0 
, 000 , 


(A24) 


where V - V2|yi| 13 cm“* AE using the parameters 
of Table S6. The effective two-level Hamiltonian in equation 
A24|is presented as equation 1 in the main text. Diagonaliza- 
tion of equation |A^ leads to delocalized eigenstates 

|P870') = cos e |P870> -I- sin 6 |P880), (A25) 

|P880') = - sin 6» |P870> -H cos 6 |P880), (A26) 

1 /2y \ 

»=-arcta„(—). (A27) 

The P870 band energy does not shift to first-order in 
but has a modified dipole strength 

|f^/>87ol^ “ I ‘'tts 6>Dp870 -H sin0Dg88oP 

= |£f 87 oI^ cos^ 0 + 2% ■ Z)b88o} COS 0 sin 0 

+ |77b8SoI^ sin^ 0 (A28) 

^ |.Dp87oP + 2^R ■ Dpsso} (^) (A29) 

\ 

(1-3 sin^ Ay) . 

(A30) 


shows that the P870 band increases its absorption intensity 
for small Ay (dipoles nearly tangent to the ring) and <I)(1) -i- 
y (dimerized unit cell). Hence, the uphill energy landscape, 
AE > 0, which sets transfer-to-trap as the bottleneck in the 
core complex, increases direct absorption in the P870 band. 
Finally, note that, in accordance with the dipole sum rule, the 
B880 band strength 

= |77b88oI^ cos^ 0-2% {Z7p870 ' -Dfisso} cos 0 sin 0 
+ |77f87oI^ sin^ 0 (A31) 

decreases by an amount such that the total absorption intensity 
I^B88ol^ l^/> 87 ol^ - (M + N)d^ remains constant. 


Equation|A29|is equation 2 in the main text, and equation A30 




b. Excitation Popuiation Redistribution 

Upon broadband excitation over the entire absorbing range, 
the initial state of the combined RC-t-LHl system is a sta¬ 
tistical mixture of the exciton states |P870'^ and |P880'^ 
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weighted by their respective dipole strengths 

+ |D^88o|2|B880')(b880'|). (A32) 


Hence, the population in the special pair is calculated, using 
equation |A25||A26| and |A28j to be 

Pp(0) = 7’r{p(0)|P)<P|) (A33) 

“ TTiTf—cos 40)+ 

4(M + N)d^[ 

+ 2% {£>p 87 o ■ £>b88o} sin 40+ 

+ I^ii88ol"(l-cos40)| (A34) 

where Tr is the trace operation. To explore the population re¬ 
distribution in general circular aggregates and P pigment clus¬ 
ters, we substitute ID^ggoP = y, |£>p 870 p = N and Ay = 0 

(which implies !R |DJ,g 7 Q ■ Z)b88o} = d^)- Expanding 

equation |A34| to second-order in 0 with this simplification 
gives immediately equation 3 in the main text. Increasing the 
number of pigments on the ring increases its radius linearly in 
M. The coupling between ring and RC states is then 

where Rq = R/32. Fig. 4 in the main text plots the percent 
increase in R-population due to extended delocalization 


APpiO) = 


(Pp(0) - N/(M + A) j 
N/(M + N) 


X 100 


(A36) 


for varying M and N using equation A34 One sees that ad¬ 


ditional E-excitation is optimized when more harvesting units 
surround fewer charge separation pigments, M > N. This can 
be understood from the second-order term in equation 3 of the 
main text. 


c. Inhomogeneous broadening 


As discussed in Supplement A3 the physiological envi¬ 
ronment of light-harvesting complexes is highly disordered 
and dynamic, leading to homogeneous and inhomogeneous 
spectral broadening mechanisms. Our aim in this section is 
to show that the simple two-level model reproduces the full 
Monte Carlo simulations reported in Supplement |A 3[ and 
gives further insight into the effects of inhomogeneous broad¬ 
ening. The analytical results for the modified P870 band 
dipole strength, P870 band energy and population redistribu¬ 



Figure 6. The full Monte Carlo simulations, detailed in Supplement 
|A 2[ and the two-level model of this section make identical predic¬ 
tions for the isolated RC and RC-LHl spectra. Both simulations take 
into account homogeneous and inhomogeneous broadening mecha¬ 
nisms. 


tion are reproduced here 


AE 


■"P870 


AE 




(A37) 


= \Dpm\^ cos^ 0 + 2% ■ 77b88o} cos 0 sin 0 


+ |7)b88oI^ sin^ 0 


(A38) 


PpiO) = 


r||77p87ol^(3 


+ cos 40)+ 


4(M + N) cP 
+ 2^R {Opg^Q ■ £>b88o} sin 40+ 

+ |7)b88oI^(1 - cos 40)1. 


(A39) 


Homogeneous broadening mechanisms can be taken into ac¬ 
count by dressing stick spectra with an appropriate lineshape 
function. The spectrum and the population of the special pair 
are averaged over the inhomogeneities of both B880 and P870 
energies 


<A(e)) = JJ^ dEpsjo dEpsm G{Eps^^)) G(£'b88o) 

X A(e, Epgjo, Ebsso) (A40) 

{Pp(0)) - JJ^ dEp^jo dEssm G{Ep^jo) G{EBm) Pp(0) 

(Ml) 


where 


A(e, £p870, Eb88o) - \D'psiQ?'g(+ - p'psioP (A42) 

is the homogeneous absorption spectrum of the P870 with a 
lineshape g(e), taken to be Gaussian in this study due to better 
agreement with our experimental data. Note that the Gaus¬ 
sian distributions, G(e), correspond to the real inhomogeneous 
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spectral widths. They are therefore different then the distribu¬ 
tions of site energies used in the full numerical simulations of 
Supplement |A 3| which are larger than the spectral widths due 
to exchange narrowing effects .El 

Figure shows how the P870 band absorption is enhanced 
when extended delocalization over LHl pigments is taken into 
account. This prediction matches that of the full numerical 
model of Supplement |A 3| Hence, this simplified model can 
be used to make further quantitative predictions. We use it 
now to understand the effect of inhomogeneous broadening 
on the excitation population redistribution. Figure shows 
that population is redistributed from antenna to P pigments for 
all realistic levels of disorder up to the full spectral linewidth 
(with associated standard deviation std ~ 200-300 cm“'), 
although the inhomogeneous broadening plays a non-trivial 
role. Even though the full expressions are used in this cal¬ 
culation, the perturbative expression, equation 3 of the main 
text, illustrates the physical principles. The sign of AE deter¬ 
mines the direction of the redistribution: AE > 0 redistributes 
population from antenna to P pigments, and vice versa for 
AE < 0. Hence, in situations where the inhomogeneous spec¬ 
tra of LHl and P overlap, there will be realizations where 
dipole strength is redistributed from antenna to P pigments 
and vice versa, perhaps leading to a decrease in the effect 
on average. However, the population redistribution is largest 
whenever AE ^ 0, so there is a non-trivial interplay between 
small inhomogeneity, which leads to a small but positive re¬ 
distribution for each realization, and large inhomogeneity, 
which leads to a possibly large (for realizations which are on 
resonance), but sometimes negative, redistribution. Interest¬ 
ingly, this interplay leads to a maximum population redistri¬ 
bution, as shown in Figure when + stdp - AEjl, 

where AE ^ 130cm ' is the energy difference between the 
P870 and B880 spectral peaks. Remarkably, the broadenings 
taken in Supplement |A 3[ which must be used to reproduce 
the experimental dark-minus-light spectra, satisfy this con¬ 
dition. Note that the inhomogeneous broadening of the RC 
taken here is equal in magnitude to that of RC models deter¬ 
mined by photon-echo and hole burning experimentJ^l This 
result gives further quantitative insight into a general design 
principle, where environment engineering may allow one to, 
as nature has, exploit this noise-assisted procesJ^. 


3. Simulating Difference Spectra: Numerical Simulations 

Illumination of the core complex induces excitation of the 
special pair. The excited state, P*, initiates a chain of charge 
transfer steps within the RC pigments which results in the oxi¬ 
dation P —> After such charge separation is produced, the 

oxidised special pair P^ no longer has the electronic proper¬ 
ties of P. In particular, it no longer possesses a transition asso¬ 
ciated with absorption at 870 nm. As a result, both the optical 
response and the coherent interactions between the P870 tran¬ 
sition and other pigments are suppressed. Additionally, the 
charge separation changes the charge distribution within the 
entire RC, which induces several effects in the optical spectra 



50 100 150 200 

stdiHi (cm''') 



stdLH! (cm (stdp fixed at 20 cm ') 


Figure 7. Averaging the change in P pigment excitation population, 
using equation |A28[ over Gaussian distributions in (A) energies of 
LHl and P exciton states and (B) energy of LHl exciton state while 
fixing the inhomogeneous linewidth of P870 at 20cm“’. Population 
is always redistributed from antenna to P pigments, but the redistri¬ 
bution depends non-trivially on the inhomogeneous linewidths. A 

maximum appears when + stdp is commensurate with half 

the energy splitting between P870 and B880. 


of the B and Hence, “light” spectra - recorded from sam¬ 
ples being subject to high illumination - will not show a con¬ 
tribution from the P870 band while “dark” spectra - obtained 
from samples subject to low illumination, such that P^ ^ P 
occurs before the absorption of the next photon - will preserve 
the optical response of P. Experimentally recorded spectra 
and the theoretical calculation for these two situations is pre¬ 
sented in Fig. 2A of the main text. 

The positions of peaks in optical absorption spectra cor¬ 
respond to the energies of optically active transitions, and 
peak broadening arises due to the interaction of excited state 
degrees-of-freedom with the environment. This interaction 
can induce decay to the electronic ground states or can re¬ 
sult in variation of the electric field in the neighbourhood of 
pigments, resulting in dynamical changes of the excited state 
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energies. Both mechanisms broaden the absorption lineshape, 
but the latter dominates in condensed matter systems. De¬ 
pending on the timescale on which variations are produced, 
the broadening is termed homogeneous or inhomogeneous. 
Homogeneous broadening acts within the lifetime of elec¬ 
tronic excitation while the inhomogeneous component leads 
to differences among elements in the ensemble at a given time. 
In this section we perform full Monte Carlo calculations to 
simulate the experimental difference spectra. This procedure 
yields very good agreement with experimentally determined 
difference spectra (Fig. 2 of the main text). 

Physiological conditions, and the associated high tempera¬ 
tures, lead to rapid fluctuations of local electric fields, mak¬ 
ing the homogeneous contribution the most relevant for our 
analysis. Accordingly, a set of parameters has been chosen to 
dress stick spectra resulting from stochastic realisations of in¬ 
homogeneities, (Tf and cTj, in the site energies and couplings, 
respectively, with a lineshape function such that the homoge¬ 
neous contribution dominates over the inhomogeneous. The 
homogeneous line shape illustrates microscopic details of the 
interaction between pigments and the surrounding protein. 
Spectra are calculated by binning the histograms of dipole 
strength from realisations of equation A5 followed by dress¬ 
ing with a suitable homogeneous lineshape function, which 
depends on the details of the transition of interest. 

The most notable differences between dark and light RC 
spectra of Fig. S4A, besides the disappearance of the P870 
resonance, are (1) slight energy shifts and (2) an increased 
height of the B and H resonances at approximately 800 and 
750 nm, re^ectively. Both changes have been previously 
documentedSl and can be taken into account theoretically by 
blue and red shifts of -H 33 cm ’ and - 39 cm ' in B and H 
energy, and a corresponding increase in the transition dipole 
moment + 0.2 D and -H 0.6 D, respectively. The electronic RC 
Hamiltonian has been modified from that of Rb. sphaeroide^ 
in order to match the experimentally determined absorption 
spectra of the RC in R. rubrum. Tables S2-S4 summarize pa¬ 
rameters for the RC used in full numerical simulations. Mea¬ 
sured and calculated absorption spectra of purified RC are 
shown in Fig. S4. 

It is well established that, under different stages of LHl as¬ 
sembly, the transition dipole strength of the B880 band, cor¬ 
responding to the Qy transition dipole moments, varies due to 
interactions with the and By transition dipole momenti^^li^. 
This increases the effective Qy transition dipole moment, and 
leads to a total absorption intensity in the core complex which 
is greater than that expected from BChl pigments in solution. 
This change of intensity is captured by an increased transi¬ 
tion dipole strength for chromophores, namely, from 6.3 D to 
^ V^6.3D = 9.8D. 

To calculate spectra for the core complex, interactions be¬ 
tween RC pigments and LHl pigments are calculated using 
the dipolar couplings given by equation A3 Fig. S4B shows 
“light” and “dark” spectra of the core complex. The theoreti¬ 
cal and experimental curves in Fig. 2 of the main text cotTe- 
spond to the subtraction of “dark” minus “light” spectra in Fig. 
S4A and S4B. Although absorption in the core complex is de¬ 
localized across RC and LHl, one can also access numerically 



A (nm) 


Figure 8. “Dark” and “light” spectra in red and blue, respectively. 
In (A) isolated RC spectra and in (B) core complex spectra. The 
difference spectra of (B) is Fig. 2(b) in the main text, also shown 
here (dark green) for clarity. The LHl complex is responsible for 
the large absorption maximum at 880 nm in (B). Theory is shown in 
continuous and experiment is shown with circles. 


the absorption spectrum of RC-like states (those states with 
more than 50% RC character), which shows a 78 % enhance¬ 
ment in absorption in the P870 region, shown in Fig. IB of the 
main text. This absorption redistribution leads to a population 
redistribution, using equation | A33 1 Taking into account delo¬ 
calization between RC and LHl pigments leads to an average 
initial excition population in the special pair PpiO) = 0.050 as 
compared to PpiO) = 0.028. Given that these numerical simu¬ 
lations account for all the RC pigments, they differ from those 
required for the initial condition of the rate model presented 
in Figure 3 since the latter model only accounts for the special 
pair in the RC. 


a. Effect of inhomogeneous broadening 

From linear spectra alone, it is difficult to discriminate con¬ 
tributions from inhomogeneous and homogeneous broaden¬ 
ing mechanisms. For the RC we chose a set of parameters 
for diagonal and non-diagonal noise, such that the standard 
deviation of the P870 resonance is about ~70 cm“'. A study 
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which combined hole burning and photon echo experiments to 
estimate the homogeneous contributiorP found, from careful 
analysis of pump-probe data, a standard deviation in the P865 
of Rb. sphaewides (analogous to P870 in R. rubrum) of 73 
cm“'. For the ring structure, 2D non-linear spectroscopy has 
provided evidence, due to the circular line-shapes observed 
in the rephasing direction of diagonal peaks at room temper¬ 
ature, of a major contribution of homogeneous broadening 
mechanisms in the LH2 aggregated. LH2 complexes have 
often been used as analogue structures of LHl for structural 
analysis. The absence of wavelength dependence in room- 
temperature fluorescence of R. rubrum also implies a domi¬ 
nant homogeneous broadening mechanisnd. However, be¬ 
yond the fact that homogeneous broadening is dominant, ex¬ 
act numbers for each of the broadening mechanisms are hard 
to read out from these experiments. Moreover, measurements 
performed at lower temperatures present major contributions 
from inhomogeneous broadening^. Here we study the robust¬ 
ness of the dipole redistribution for different degrees of inho- 
mogenoeus noise in order to clarify their non-trivial effects. 

Figure 1^ presents the enhancement of the P870 resonance 
due to the coherent coupling with ring pigments for different 
levels of inhomogeneities in the intra-ring couplings, i.e. the 
non-diagonal Hamiltonian elements, while keeping the diag¬ 
onal disorder fixed. Off-diagonal disorder has been shown to 
be crucial to explain fluorescence anisotropjHl. For compari¬ 
son, this Figure also shows the isolated RC spectra obtained 
experimentally. In each case the isolated RC spectrum cal¬ 
culated numerically is fitted to the experiment by tuning the 
homogeneous broadening. This study reinforces the findings 
of Supplement A 2 c That is, the dipole strength redistribu¬ 
tion is qualitatively robust to inhomogeneous broadening, and 
the redistribution is maximized for intermediate levels of dis¬ 
order. 


4. Generalization to other geometries in purple bacteria 
species 

The closed-ring circular symmetry of R. rubrum seems to 
be robust in that specieJ^, however, different geometries have 
been observed in other species. For instance, there is good 
evidence, from high-resolution atomic force spectroscopy, as 
well as cryoelectron spectroscopy, that the LHl ring of the 
PSU from Rb. sphaeroides is interrupted by other protein 
subunits, corroborated by biochemical evidence that only 24 
BChls are present in the LHl structur^i^ More recently, the 
structure of T. tepidum, the highest-resolution structure yet ob¬ 
tained, has provided evidence of an oval-shaped LHl with an 
eccentricity of 0.4^^. We have performed detailed numerical 
simulations for structural perturbations such as these in order 
to prove the robustness of the RC cross-section enhancement 
in situ, as compared to the isolated forms. 

We perform Monte Carlo simulations of Rb. sphaeroides 
with parameters based on pr eviou s agreement between sim¬ 
ulated and measured spectr^SEH in our analysis we use 
the structure of Ref.ES, the LHl parameters from Ref.® and 
tune the site energies to obtain an LHl absorption spectra 


which peaks at 875 nm, as observed experimentally. The IR 
resonances in Rb. sphaeroides have been measured to vary 
betwee n 87 0 and 875 nm, depending on the reconstitution 
proces ^221211 \Ye develop calculations based on the ring res¬ 
onance peaking at 875 nrrP^. Hamiltonian and noise param- 
eters from the RC of Rb. sphaeroides are taken frorrpSi while 
RC homogeneous broadening parameters are presented in Ta¬ 
ble S5. As a check on parameters, spectra of RC and LHl for 
this species are calculated and presented in Fig. they are 
qualitatively similar to previous worlsPH, and most imporantly, 
the anti-symmetric P band and the LHl band peak around 865 
and 875 nm, respectively. 


In a first set of simulations, the complete 16a/3 circular LHl 
aggregate is investigated. Figure [T0|A) shows the enhance¬ 
ment in the absorption strength of the P band due to the co¬ 
herent interaction with the LHl pigments. The enhancement 
is less than that observed in R. rubrum because, here, the LHl 
dipole moment is smaller. This leads to a reduced coherent 


coupling and dipole redistribution according to equation A30 


As noted earlier, hyperchromism in R. rubrum leads to inten¬ 
sity borrowing of the Qy bancf^. Such an effect has not yet 
been studied in any other species, and hence cannot be gen¬ 
eralised or disregarded. Accordingly, an exact value for the 
dipole strength associated with the Qy transition dipole of pig¬ 
ments in other species is hard to obtain; values of 6.1, 7.7 and 
11 D have been reported^^I^. Hence, for these simulations 
6.3 D (dashed), corresponding to the bare transition dipole of 
BChla, and 9.8 D (dotted), corresponding to the measured hy¬ 
perchromism of R.rubrum, are used in Figure [T0| A). With 6.3 
D the increase of the RC cross-section is 17 %, while 9.8 D 
yields a 75 % increase. 


The second set of simulations takes into account more re¬ 
cent structural information that displays an LHl complex ar¬ 
ranged on an arc with 14 afi units, corresponding to 28 pig¬ 
ments. For comparison purposes, the parameters of the 16- 
fold symmetric case are used. The absorption spectra of the 
LHl complex for the 16-fold symmetric and 14-dimer rings 
are presented in Figure [T0|B). The missing dimer structure, 
shown in the inset, produces a spectrum which is slightly 
wider than that of the complete circular aggregate. Hence, 
in order to compare changes in the RC spectrum we reduce 
the inhomogeneous broadening to cry = 30 cm“^ in the 14- 
dimer case such that its corresponding spectrum is as wide as 
that from the 16-fold symmetric structure. In addition, no hy¬ 
perchromism is posited, so that the bare dipole moment, d = 
6.3 D, is used. Figure [T0|C) shows the enhancement in the 
absorption strength of the P band due to the coherent inter¬ 
action with the LHl pigments. Note that the enhancement 
in the 16-fold symmetric case is greater than that of the 14- 
dimer case. As shown in equation A30 the magnitude of the 
dipole strength enhancement depends linearly on the number 
of dipoles in the special pair and LHl ring, N and M, respec¬ 
tively. It is therefore natural to expect a smaller enhancement 
with N - 2, M - 28 as opposed to N - 2, M - 32. To avoid 
this pigment number effect and investigate the effect of mi¬ 
croscopic details of the structure. Figure [TotP) compares RC- 
like spectra of the 16-fold symmetric ring and the 14-dimer 
ring when the transition dipole of the latter is increased from 
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Figure 9. The effect of inhomogeneous broadening in the P870 band. In all plots experimental “dark” spectra for isolated RC are represented 
with circles; they are almost indistinguishable from numerical calculations for isolated RC (solid blue lines). In black are shown the spectra 
of RC-like states, i.e. states with an excitation population of more than 50% in the RC. In A, B and C the standard deviations of both 
nearest-neighbour ring couplings and P pigments energies, are 20, 100 and 200 cm“', respectively. These figures show that the dipole strength 
redistribution is qualitatively robust to noise and maximized for intermediate levels of disorder. 


6.3 D to 6.3 X 32/28 D. As apparent from Figure [T0|^D), this 
consideration explains most of the difference between the two 
spectra. In order to understand why the curves do not overlap, 
further comment is required. 

In general, perturbations from a perfect ring lead to a re¬ 
distribution of transition dipoles in the ring, to other levels 
different than k - ±1. In the case of the missing dimers, the 
k - 0 level gains dipole strength at the expense of a concom- 
mitant reduction in the k - +l band. Given that k - 0 level 
is further off-resonance than the k - ±l from the P870 band, 
its dipole strength does not appreciably increase the RC cross 
section. 

A preliminary investigation of oval perturbations shows 
similar results to Figure[T^ Oval-like structures have been re¬ 
ported f or the PSUs of the purple bacteria T. tepidum and Rps. 
palustri^^^. We have observed that the quantitative enhance¬ 
ment can be closely related to the dipole strength distribution 
in the LHl manifold, itself dependent on the geometrical con¬ 
figuration. 

This study shows that a proper quantification of the dipole 
strength redistribution, from antenna to P pigments, can give 
insight into the fine details of the couplings, energy landscape 
and ring geometry. However, the most important finding of 
these examples is that we always find an increase in the RC 
cross-section due to the coherent interaction with the symmet¬ 
ric or the slightly-perturbed ring, which supports the assertion 
of generality across purple bacteria core complexes. 

5. Methods for determination of the RC-P870 contribution 
after oxidation of the intact, membrane embedded PSD with 
stoichiometric amounts of K 3 Fe(CN )6 

a. Diagnostic features of the RC redox spectrum 

It is now well-accepted that the oxidation of the purified 
RC with small amounts of a strong oxidizing agent, such 
as K 3 Fe(CN) 6 , leads to spectral changes in the near-IR re¬ 
gion which correspond perfectly to those observed after pho- 
tobleaching with actinic light. A typical example, obtained 
by the K 3 Fe(CN) 6 -induced oxidation of the purified RC (in 


detergent solution) from R. rubrum is shown in Fig. S7A. Un¬ 
der these conditions of oxidation the P870 band is essentially 
abolished, and the 803 nm peak, corresponding to the acces¬ 
sory BChl is slightly increased in intensity and blue-shifted. 
It has been shown previousl}® that the absorption transitions 
contributing to P870 are in fact split upon oxidation to P870 h-, 
with one component being shifted to about 800 nm (hence the 
observed changes to the 803 nm peak) and the other to 1260 
nm (which is not visible with our instrumentation). For di¬ 
agnostic purposes with real biological preparations, the redox 
difference (reduced RC - oxidized RC) spectrum (Fig. S7B) 
is a more reliable indicator of RC activity, since any contribu¬ 
tions due to inactive RC are cancelled out. 

A very useful feature of the RC redox spectrum is the 
crossover point and relative intensities of the flanking regions 
(indicated by a dotted line in Fig. S7B) at the 803 nm region. 
This feature is uniquely due to RC oxidation, and is not ob¬ 
served from isolated LHls subject to oxidation. In a complex 
mixture, therefore, this feature can be used as a diagnostic 
feature to distinguish RC and LHl contributions, respectively. 


b. Measurement conditions for obtaining spectra of 
purified RCs and chromatophores 

Chromatophores (intracytoplasmic membranes containing 
the PSUs and reaction centers) were isolated from anaero- 
bic, phototrophic cultures of R. rubrum S1 as describecf^E^l^. 
Absorbance spectra of 500 ^d samples (chromatophores; di¬ 
luted in 20 mM TrisHCl, pH 8.0 (20T8)) were usually mea¬ 
sured at room temperature in 2 mm path-length cuvettes with 
a Jasco V-560 UV/VIS spectrophotometer equipped with a 
near-IR-sensitive photodiode detector for turbid samples (slit 
width 2 nm, scanning speed 200 nm/min, photomultiplier re¬ 
sponse fast). Purified RCs (purifed by a slightly modified 
method of Vadebonceur et al™, were diluted in 20 mM Tr¬ 
isHCl, pH 7.7 containing 50 mM NaCl and 0.05 % (w/v) 
dodecyl-jS-D-maltoside, prior to measurement. For routine 
chemical oxidation of isolated RCs (A870 (2 mm path-length) 
= 0.016), a few grains of potassium ferricyanide (K 3 Fe(CN) 6 ) 
were added directly to the cuvette prior to measurement. 
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Figure 10. The effect of breaking the 16-fold ring symmetry on 
the dipole redistribution. Figure (A) shows the RC (continuous) 
and RC-like for hyperchromic (d = 9.8 D) (dot-dashed) and non- 
hyperchromic (d = 6.3 D) (dashed) spectra in the 16-fold sym¬ 
metric LHl structure. Figure (B) compares the LHl-like spectra 
(spectrum of states with more than 50% excitation population in 
the LFIl pigments) for the 16-fold symmetric (continuous) and 14- 
dimer (dashed) structures using the same amount of inhomogeneous 
noise. Over the 16-fold symmetric curve lies a 14-dimer absorption 
spectra (dot-dashed) with reduced inhomogeneous width. In (C) a 
comparison between RC-like states of the 16-fold symmetric struc¬ 
ture (dashed) and 14-dimer structure with reduced inhomogeneity 
(dotted-dashed). In (D) the same comparison is shown, but the dipole 
strength in the LHl for the 14-dimer case is 7.2 D instead of 6.3 D to 
account for the effect of the missing pigments. Most importantly, one 
sees an appreciable increase in the dipole strength of the P pigments 
when extended delocalization is taken into account. 


c. Determination of the RC concentration from a spectrum 
obtained from totai chromatophores 

A precise determination of the RC concentrations is an es¬ 
sential prerequisite for stoichiometric oxidation. In princi- 



Wavelength (nm) 



Figure 11. (A) Absorption spectra of reduced and oxidized RC, re¬ 
spectively. The position of the special pair (P870, Gaussian fit: 

= 869 nm, cr = 19.6 nm) is indicated. (B) Redox (reduced RC - oxi¬ 
dized RC) difference spectrum. The position and intensity of the 803 
nm crossover, used as a reference point for fitting, is indicated. 


pie, the RC concentration in total R. rubrum chromatophores 
might be determined from the 803 nm peak, which is ex¬ 
clusively due to the RC-accessory BChl. Unfortunately, the 
881 nm peak of the LHl exhibits significant intensity in the 
flanking region at 803 nm (see Fig. S8), which makes the 
assignment of the baseline solely due to the RC 803 nm 
peak difficult. To solve this problem unambigously, we com¬ 
pared the spectrum from the wild-type (LHl-i-RC) with that 
of a RC-minus mutant, SKALM, where the L- and M sub¬ 
units of the RC have been deleted (described in Ref.l^. We 
have previously shown that the absorption spectra and also ul¬ 
trafast transfer times for energy transfer amongst BChl and 
carotenoid pigments of the LHl of the mutant are unchanged 
when compared with those of the wild-type LHlISI. We then 
fitted the SKALM spectrum to the wild-type spectrum, using 
the criterium that the right flanking region, between 882 - 900 
nm, is due solely to the contribution from the LHl. This fit 
(see Fig. S8, expanded insert) enabled us to set a reference 
point for the position of the RC baseline. We note that the 
SKALM fitted spectrum now corresponded perfectly to the 
position of the LHl band at 589 nm (data not shown). This 
calculation showed that for the LHl-i-RC chromatophores, 
about 48% of the intensity at 803 nm arises from the RC ac¬ 
cessory BChl. Using the millimolar extinction coefficient of 
e(803 nm) = 329 mM^'.cm^' we calculated that 243 nmol 
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Figure 12. Absorption spectra used for the determination of the RC 
content. The effective RC baseline was determined from the differ¬ 
ence in intensity at 803 nm from spectra of isolated chromatophores 
from the wild-type SI (LHl -l- RC, brown line) and a RC- mutant 
SKALM (LFll only, green line). The calculated RC spectrum (us¬ 
ing the spectrum obtained from purified RC (see Fig. S7A) is shown 
in blue. The spectrum shown in red corresponds to the (LHl-l-RC) 
chromatophores minus the calculated, purified RC spectrum. 


RC are present in the cuvette under the conditions used. 


d. Stoichiometric oxidation of the RC with KiFe(CN)(, 

For a precise titration of the RC special pair redox state 
in chromatophores (A881 (2 mm path-length) = 0.37), in¬ 
cremental amounts of K 3 Fe(CN )6 (0 - 560 nmol (stock so¬ 
lution 20 mM K 3 Fe(CN )6 solution (in H 2 O))) were added to 
a chromatophore suspension, containing approximately 243 
nmol RC, as calculated from the extinction at 803 nm. Gener¬ 
ally, the chromatophores were incubated with redox agent at 
25°C for 5 min prior to measurement in the Jasco spectropho¬ 
tometer, as described above. Prior to quantitative analysis, the 
spectra were normalized to the baseline intensity at 640 nm, 
where no absorption due to pigment occurs. 

The results of a typical redox titration are shown in Fig. 
S9A. We observed, that even at stoichiometric levels of 
KjFe(CN)(, to RC, that the 882 nm peak due to the LHl Qy 
absorption maximum is also reduced (note that this bleaching 
corresponds to chemical oxidation of the LHl BChl), and this 
reduction is concomitant with that of the special pair oxida¬ 
tion. To our knowledge, this is the first time that the effects 
of stoichiometric oxidation upon the LHl complex have been 
reported. 

To proceed further, we made the initial assumption, that at 
higher, but still near stoichiometric quantities of K 3 Fe(CN) 6 , 
that the oxidation of the LHl BChls (measured as the reduc¬ 
tion of peak height at 881 nm) should be linear, whereas that 
of the RC (measured at 803 nm) should be constant. This can 
be most easily visualized by calculating the difference spec- 



Figure 13. Reduction of the near-IR (LHl-t-RC) peak intensity fol¬ 
lowing stoichiometric additions of KsFefCNle. (A) The primary 
spectral data, obtained from chromatophores equilibrated with dif¬ 
ferent concentrations of KsFefCNle (given in nmoles). The initial 
data set was normalized to the intensity at 640 nm, where no ab¬ 
sorption due to pigment occurs. Inset: enlarged near-IR region of 
the spectrum. (B) First near-IR difference spectra (AAbs), obtained 
from the data shown in (A). Each spectrum corresponds to the ditfer- 
ence between the spectrum obtained at a given amount of KsFefCNle 
(indicated) minus the reference spectrum (taken in the absence of 
K 3 Fe(CN) 6 ). Inset: the spectral difference obtained over the com¬ 
plete range. Note that the absorption change due to KsFetCN)^ is 
easily visible (arrow), and can be used for control purposes. 


trum (A881 (0 mM K 3 Fe(CN) 6 ) - A881 (x mM K 3 Fe(CN) 6 ) = 
(AAbs881)). The spectral results, shown in Fig S9B, confirm 
our hypothesis: above about 240 nmol K 3 Fe(CN) 6 , no further 
change in the 803 nm peak is observed with increasing oxi¬ 
dizing agent, but the A881 nm peak is incrementally reduced 
further. The inset in Fig. S9B also shows that no spurious 
effects to turbidity changes arise during this process. Quan¬ 
titative analysis of these data (Fig. SlOA and B) confirm the 
linearity of LHl oxidation at K 3 Fe(CN )6 levels where the RC 
is completely oxidized (Fig. SlOA). The usable range for fur¬ 
ther analysis, where the RC oxidation is constant (Fig. SlOB) 
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Figure 14. Variation of the peak amplitudes of the peak intensities 
at 881 nm (A) and 803 nm (B) respectively, obtained from the dif¬ 
ference spectra shown in Fig. S9B. In (A) the dashed line was de¬ 
termined by linear regression, using the last four points (320 - 560 
nmol) as input data. Flowever, inclusion of the data point at 240 
nmol had almost no effect upon the regression curve. Note that the 
intensities shown have been re-scaled for presentation clarity. 


was shown to be between 400 - 560 nmol K 3 Fe(CN) 6 . 

Thus, the LHl spectral contribution present in the spectra 
obtained form oxidized LHl-i-RC chromatophores should be 
obtainable by subtracting pairs of difference spectra (AAbs) 
obtained for K 3 Fe(CN )6 amounts between 400-560 nmol. In 
principle, any of the pairs might be employed for this pur¬ 
pose, but we have found it convenient to obtain a series of 
difference-difference spectra by subtracting sequental pairs 
(i.e. AAbs320-AAbs400, AAbs400-AAbs480, AAbs480- 
AAbs560). With this series in hand (Fig. SllA), we need 
to perform one final control to prove that each member of this 
series corresponds to the same spectral form. This is done 
by a simple transformation. First, each (now noisy) AAAbs 
spectrum is normalized to the same maximal intensity (Fig. 
S1 IB). If the normalized spectra correspond to the same spec¬ 


tral form, then a plot of AAAbs values for each pair should fall 
upon the same straight line. Fig. SI 1C shows that this expec¬ 
tation is fulfilled. Therefore, we can now employ the parame¬ 
ters of the derived LHl spectral difference for determining the 
contribution in the redox spectrum obtained from the lowest 
amounts of K 3 Fe(CN) 6 , where the contribution from RC and 
LHl, respectively are comparable. 

The final deconvolution was performed as follows. We 
chose to use the primary difference spectrum AAbs(160) for 
the analysis, since visual inspection suggested significant 
spectral contributions from both LHl and RC to be present. 
The procedure was as follows. Initially, the secondary differ¬ 
ence spectrum AAAbs(560-480 nmol) (see Fig. SllA) was 
fitted to the spectrum AAbs(160), using the constraint that the 
least-squares difference at the far red flank (which arises al¬ 
most exclusively due to the LHl component) must be min¬ 
imized. The result, AAbs(160) - AAAbs(560-480 nmol), is 
shown in green in Fig. S12A and B, and the fitted LHl com¬ 
ponent is shown as a dashed cyan line in Fig. S12A. Also 
shown is a Gaussian fit {Amax - 868 nm, cr= 19.6 nm) to the 
“P870” peak of the final spectrum, which is thus 1 nm blue- 
shifted compared to that obtained for the purified RC in so¬ 
lution (shown as the orange dashed line in Fig. S12B). In 
Fig. S12A, another secondary difference spectral component 
(AAAbs(320-240 nmol) is also shown (brown dashed line)). 
The same fitting as above, repeated using this latter spectral 
component, led to an identical result (not shown). In the main 
text, the spectrum shown in Fig. 2A has been calculated using 
the difference AAbs(160)-AAAbs(320-160 nmol). The final 
result however, is essentially identical to that shown in Fig. 
S12A,B. 


e. Methods for core complex photobleaching 

The measurement of RC photobleaching was performed 
by comparing the “dark” spectrum obtained from total chro¬ 
matophores (A881 (2 mm path-length) = 1.06) to that ob¬ 
tained during illumination with actinic light (using a 150 W 
halogen lamp). The actinic light was filtered through 3 cm of a 
1 M CUSO 4 solution. The amount of light quanta reaching the 
cuvette during illumination in this experiment was 40 nmol 
crn^^.s^^ (or 14000 Lux). In the dark-light experiments, the 
internal photomultiplier of the Jasco spectrophotometer was 
used instead of the photodiode, because it could be protected 
more easily from any stray actinic light by a long-pass filter 
(Asahi, XVL 0610, 1.5 mm) and a self-built chamber of black 
cardboard. To improve the signal-to-noise-ratio, the scanning 
speed was reduced to 100 nm/min, and the response was set to 
slow. In this experiment, the contribution to LHl fluorescence 
was negligible. This was confirmed by performing control 
experiments using an Avantes diode array spectrophotometer, 
where the intensity of the measuring beam was sufficient to 
induce measurable LHl fluorescence. 
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Figure 15. (A) The secondary difference spectra, obtained from the 
data shown in Fig. S9B. Each spectrum corresponds to the difference 
taken from two sequential pairs e.g. AAAbs(A560-480) = AAbs(560) 
- AAbs(480). (B) The difference spectra in (A) normalized to the 
same intensity value at 881 nm. (C) Correlation plot of the AAAbs 
intensities from sequential spectral pairs (e.g. AAAbs(560-480) vi. 
AAAbs(480-400)). The intensities shown have been re-scaled (x 10“) 
for presentation clarity. Only the intensity range corresponding to the 
wavelength range 800 - 900 nm has been plotted. 


6 . Charge Separation Dynamics 

Excitonic transfer between the RC and LHl occurs via a 
Coulomb exchange mechanism on a picosecond time-scald^, 
while vibrational dephasing destroys coherences in hundreds 
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Figure 16. Deconvolution of the primary difference (AAbs) spectra 
into LHl and RC components. (A) The primary difference spectrum 
AAbs(160) (blue line), obtained at 160 nmol KsFelCN)^, which con¬ 
tains about equal intensities from both LHl and RC, was used for 
analysis. The fitting procedure was performed using the secondary 
difference spectrum AAAbs(560-480 nmol) (see Fig. SI lA) (see sec¬ 
tion SI (E) for details). The result, AAbs(160) - AAAbs(560-480 
nmol), is shown in green in (A) and (B), the fitted LHl component 
(dashed cyan line in (A)), and a Gaussian fit = 868 nm, cr = 
19.6 nm, dark blue dashed line) to the “P870” peak of the final spec¬ 
trum, which is thus 1 nm blue-shifted compared to that obtained for 
the purified RC in solution (shown as the orange dashed line in (B)), 
are also shown. In (A), another secondary difference spectral compo¬ 
nent (AAAbs(320-240 nmol)) is also shown (brown dashed line). The 
same fitting as above, repeated using another secondary difference 
spectral component (AAAbs(320-240 nmol), brown dashed line) led 
to an identical result (not shown). Note that the intensities shown 
have been re-scaled (x 10^) for presentation clarity. 


of femtosecondJ^mi The Coulomb interaction de-excites an 
initially excited electron in the donor complex while simul¬ 
taneously exciting an electron in the acceptor complex. As 
dephasing occurs, the excited states of donor and acceptor 
acquire random relative phases, which inhibits quantum me¬ 
chanical delocalisation. 

Transfer rates measured from pump-probe experiments 
agree with those obtained from calculations using general¬ 
ized Forster theorjG^. The theoretical approach to date aver¬ 
ages transfer from different excitonic delocalized states from 
a donor complex, assumed to be in a canonical equilibrium 
of excited states, to the acceptor complex. LHl—> RC trans- 
















pp{Q) = Pp{Q) 
Plwi(O) = 1 - p/>(0) 
Pp+(0) = 0 
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Figure 17. Increase in efficiency due to quantum delocalization as a 
function of the relevant incoherent tranfer rates, / and /, according to 
equation |A46| Since transfer between light-harvesting units occurs 
much more quickly than excitation loss (/// ~ 6), the increase in 
efficiency is small for R. rubrum. However, this plot illustrates that 
this design principle could lead to higher increases in efficiency in 
artificial systems like quantum dots and, in general, when the system 
is limited by loss (/// < 1). 


fer occurs through the second and third lowest exciton lying 
states, giving rates in agreement w ith an experimentally mea¬ 
sured transfer time of / = 40 Back transfer from 

RC to LHl occurs in a calculated time of b = 8.1 pP, close 
to the experimentally measured 7-9 ps, which was estimated 
from decay kinetics after RC excitatiorP^l and further verified 
numerically by careful modelling of environmental effectP. 
The first electron transfer step, P* P^, occurs in the RC 
within c = 3 pP. Fluorescence, inter-system crossing, inter¬ 
nal conversion and further dissipative mechanisms, have been 
included in kinetic models with a single lifetime of I = 250 
pP. The rates are presented schematically in Fig. 3 of the 
main text. 

Excitonic transfer among LHl and RC pigments occurs in 
a regime where coherent delocalisation across RC and LHl 
shows no appreciable effectP. Therefore, to describe exci¬ 
tation transfer in the core complex a model of classical rate 
equations taking into account LHl, P, a charge separation 
state P^ and a loss sink is sufficient. The system evolves ac¬ 
cording to 


_ 1 11 

pLHi - -T Plh\ + t Pp —T Fi/ri 

/ b I 

1 111 

pp^ Plhi - tPp- -Pp- jPp 

J o c i ^^43^ 

pp* ^ - Pp 
c 

pioss = y (Plh\ + Pp)- 


The initial condition 


(A44) 


for the coherent case, PpiQ) follows from equation [A33[ while 
Pp(0) = NUN + M) when the interaction between F-pigments 
and the LHl is neglected. These equations can be solved ana¬ 
lytically for the probability of charge separation at long times: 


ppAt 


oa) - bl 


I / + fPp(0) \ 

[cfl + b(c + l)(f + l)j' 


(A45) 


We are interested in the increase of charge separation due 
to extended delocalization, which changes the initial condi¬ 
tion Pp{Q). Defining 6 to be the initial population in the P- 
pigments without extended delocalization, and 6 > 6 with 
extended delocalization, the percent increase in charge sep¬ 
aration, which is the change in efficiency Ae, is given by: 




l + f6 


(A46) 


As shown in the main text and in the SI, d = ^ and 6 ^ 
Ligure [T7| shows the increase in charge separation according 
to Equation |A46| for variable values of 1//, the ratio of the dis¬ 
sipation to the forward transfer time scales, which determines 
the increase in efficiency. Since / > / in R. rubrum, charge 
separation is very efficient (~ 83%) and the benefit from pop¬ 
ulation redistribution due to extended delocalization is ~ 1%. 
However, for lossy systems f < I, the benefit of extended 
delocalization can be very large. Ligure [T^ gives a concrete 
example of a system where the probability of charge separa¬ 
tion increases from 0.21 to 0.33, an increase in efficiency of 
54%. 


7. Higher Plants 


Photosystem 1 (PSl) and photosystem 2 ^S2) are the main 
light-harvesting photounits in higher plantPD. To explore the 
effects of extended delocalization in these systems, we exam¬ 
ine Hamiltonians of the form 


H = 


(Hrc 

I ^ 



(A47) 


where Hpc is the Hamiltonian of the full RC alone and Hpu 
is the Hamiltonian of the antenna complexes. The coupling 
between RC and LH pigments is captured by the matrix V. 
In this simplified calculation, site energies are taken from the 
literature and couplings between pigments are calculated in 
the point-dipole approximation using d - 4.48 D, /c = 1 and 
d - 4.4 D, A- = 1.5 for PSl and PS2, respectiveljISMUl Dipole 
moments, Di, are calculated along the direction connecting the 
Nb and No atoms of the porphyrin rings of each chlorophyll 
using the atomic coordinates taken from the X-ray structures 
(PDB accession codes IJB O and 3WU2), and rotated 15" to¬ 
wards Nc in the case of ps Lluctuations in the pigment 
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Figure 18. Optimizing charge separation in an artificial system. 

Tuning the resonance condition, AE = 9 cm“’, maximizes the popu- 
lation redistribution (Pp(0) = Iz = 0.29) according to equation 


A3 3 


leading to a 400 % increase in direct excitation of the P-pigments. 
For a situation in which the pigments are subject to fast losses ( f = 
40 ps, b = 8 ps, c = 3 ps, 1 = 15 ps), the population redistribution 
increases the probability of charge separation from 0.21 to 0.33, a 
54 % increase in efficiency. 


energies, drawn from a Gaussian distribution, capture spectral 
inhomogeneities, which are known to fluctuate with tempera¬ 
ture. We take Gaussian standard deviations of cr = 80 cm ' (fit 
from low temperature spectra) and cr = 28 cm“* (fit from r oom 
temperature spectra) for PS 1 and PS2, respectiveljHSES For 
each realization of the Hamiltonian, H, we And a unitary U 
which diagonalizes it and calculate the dipole moments of the 

exciton states, Ik) = 2; 10^ according to D'/^ - 2,- UikOj. 

The population in the P pigments is calculated according to; 


P(0) = 


ZJDVlaXal 

Za\D'a\^ 


N 

PF(0) = ^7’r{p(0)|P,)<P,|). 
1=1 


(A48) 

(A49) 


where N — 2 for PS 1 and PS2 monomer and N - 4 for the 
PS2 dimer. We compare the calculation of equation | A49| for 
y = 0 (uncoupled RC and antennae) and V 4 0 (coherently 
coupled RC and antennae). As summarized in figure [T^ these 
calculations predict an increase in the excitation of P pigments 
due to excitonic delocalisation over antenna chlorophylls ob¬ 
served in the X-ray structures. Excitation redistribution from 
antenna to P pigments is robust to, albeit reduced by, disor¬ 
der. A more sophisticated calculation, which we leave for a 
future work, will likely require the determination of nearest- 
neighbor couplings beyond the dipole-dipole approximation, 
which is known to break down at distances of approximately 
10 A. 
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Figure 19. Excitation redistribution from antenna to P pigments in 
PSl (A) and PS2 (B) is robust to Gaussian disorder in the site ener¬ 
gies. The average change of P population in PSl shifts from +23% 
with no disorder to -1-14% with 80 cm“' inhomogeneous broadening; 
the distribution has a Lorentzian FWHM of 20%. In the PS2 dimer, 
the redistribution shifts from +16% with no disorder to +12%, Gaus¬ 
sian std 11%, with 28 cm“' inhomogeneous broadening, and the PS2 
monomer shifts from -11% to +12%, std +15%, with the same in¬ 
homogeneous broadening. Hence, although disorder levels the dif¬ 
ference between the PS2 dimer and monomer in average, extended 
couplings in the dimer increase robustness to site energy fluctuations, 
as the distribution is narrower for the dimer. 


8 . Does quantum delocalization across photosynthetic 
pigments of the RC-LH1 complex contribute to evolutionary 
fitness? 

At first glance, it would seem that an evolutionary advan¬ 
tage due to quantum mechanical delocalization is unlikely, as 
the improvement (~1%) is small compared to the efficiency. 
Experi mentaly -derived rates yield overall efficiencies of about 
85-95‘%p2E2|43| Under our normal laboratory growth condi¬ 
tions, we employ light intensities of about 24 yuE.m^^.s^', 
equivalent to 1.5 x lO'® photons.m^^.s '. However, the com¬ 
fortable laboratory conditions do not resemble the natural eco¬ 
logical niche of P. rubrum, and only in the latter has evolu¬ 
tionary pressure molded the presently observable (LHl+RC) 
PSU structural organization. R. rubrum is a fresh-water bac- 
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terium, which is found in lakes under conditions of low light 
intensity and low oxygen. We can estimate the probable wa¬ 
ter depth of the ecological niche by reference to the oxygen 
tension which allows the expression of photosynthetic mem¬ 
branes. In a study by one of us (R. GhoslJ^ we have shown 
that only under a p02 of 0.3% are photosynthetic membranes 
expressed in R. mbrum. In a meromictic lake (where a chem¬ 
ical gradient between oxygen from the surface, and nutrients 
from the lake sediment occurs) such as Lago di Cadagno in 
the Swiss alps, a p02 of 0.3% occurs at a depth of between 
6-11 rrPl depending on the season of the year. At this wa¬ 
ter depth, the light intensity may only be 1 - 0.1 pE.m^^.s '. 
Under these conditions, taking geometrically-determined self¬ 
shading due to the clustering of the intracellular photosyn¬ 
thetic membranes within the cell, the R. mbrum PSU must 
be adapted to extremely low light intensities - the incident 
light intensities could be as low as 60 photons.nm^^.ms ', 
which may be severely limiting, considering that the time 
for one turnover of photoinduced cyclic electron transport 
is about 1 m^. Based on these considerations, it is not 
unreasonable to speculate that the increase of the effective 
extinction coefficient of the RC due to quantum delocaliza¬ 
tion is related to the evolutionary pressure determined by the 
growth light intensity conditions. Finally, note that the pre¬ 
dicted enhancement of P-pigment excitation in the species 
we studied {R. rubrum, APp(Q) - 100%; Rb. sphaeroides, 
APp{Q) ~ 20 - 75%, depending on the hyperchromism; PSl 
and PS2, APp{0) = 12 - 14%) is correlated with the light in¬ 
tensities that these organisms experience under physiological 
conditions. More explictly, lower light intensities correlate 
with a greater enhancement, supporting our claim that this ef¬ 
fect could be both general and biologically relevant. 

Appendix B: Supplementary Tables 


Table I. Parameters for the LHl of R. rubrum which capture the ab¬ 
sorption spectrum of Fig. S4. 


Parameter 

Value 

Description 

7 

10.06° 

intra-dimer angle 

A7i,A72 

19.9°, 17.6° 

dipole tilt away from tangent 

01 , <^2 

6.24°, 4.81° 

dipole tilt out of membrane 
plane 

ri,r2 

46.51 A, 47.27 A 

radius of a, p rings 

z 

0.63 A 

vertical displacement among a 
and p chromophores 

dtHi 

9.8 D 

dipole moment for Qy transition 
of LHl 

K 

2 

relative permittivity 

fff, 6/3 

12646,12656 cm^* 

chromophore energies of a and 
P site 

2 i, Qi 

600, 377 cm-' 

intra- and inter-dimer coupling 

o-j 

60 cm“‘ 

magnitude of nearest neigh¬ 
bours coupling strength 

disorder 

CTlHI 

80 cm-* 

magnitude of LHl pigments en¬ 
ergy disorder 

^LHl 

36.5 nm 

half width half maximum ho¬ 
mogeneous dressing L(/l, 4°, F) 
for LHl transitions 


Table II. Parameters for the active RC of R. rubrum which capture 
the absorption spectrum of Fig. S4. 


Parameter 

Value 

Description 

dp 

6.3 D 

dipole moment for Qy transition 
of P BChls 

dAp 

6.7 D 

dipole moment for Qy transition 
of B 

dp 

8.2 D 

dipole moment for Qy transition 
of H 

a-p 

70 cm“‘ 

magnitude of P pigments en¬ 
ergy disorder 

(Tap 

40 cm“‘ 

magnitude of B energy disorder 

cr Bphy 
c Vp 

70 cm“‘ 
50 cm“* 

magnitude of H energy disorder 
magnitude of P pigments cou¬ 
pling disorder 

Tp 

21 nm 

standard deviation homo¬ 
geneous dressing for P870 
transition 

Fa/> 

9.3 nm 

standard deviation homoge¬ 
neous dressing for B transition 
at 802 nm 

^ Bphy 

35 nm 

half width half maximum ho¬ 
mogeneous dressing for H 


Table III. Hamiltonian (in cm * units ) for the RC of R. rubrum in the 
basis {|P1), \P2 ), \APl ), \AP2 ), \BPhl ), \BPh2)]. Missing elements 
are calculated using equation | A3 1 


('11990 

500 

-60 

-70 


500 

11990 

-70 

-60 


-60 

-70 

12537 


140 

-70 

-60 


12537 

130 



140 


13 220 




130 

13 220; 
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Table IV. Parameters for the oxidised RC of R. rubrum which capture 
the absorption spectrum of Fig. S4. 


Parameter 

Value 

Description 

dp 

OD 

dipole moment for Qy transition 
ofPBChls 

dAp 

6.9 D 

dipole moment for Qy transition 
ofB 

dp 

8.8 D 

dipole moment for Qy transition 
ofBPh 

£ap 

12 537 cm-* 

B energy 


13 220 cm-‘ 

BPh energy 


Table V. Parameters for the RC of Rb. sphaeroides 


Parameter 

Value 

Description 

dp 

6.3 D 

dipole moment for Qy transition 
of P BChls 

dAp 

6.3 D 

dipole moment for Qy transition 
of B 

dp 

6.3 D 

dipole moment for Qy transition 
of H 

a-p 

240 cm“‘ 

magnitude of P pigments en¬ 
ergy disorder 

O'AP 

90 cm“‘ 

magnitude of B energy disorder 

O' Bphy 

120 cm-* 

magnitude of H energy disorder 

O' Vp 

0 cm“‘ 

magnitude of P pigments cou¬ 
pling disorder 

r,> 

16 nm 

standard deviation homoge¬ 
neous dressing G(/i,/l®,r) for 
P870 transition 

Tap 

11 nm 

standard deviation homoge¬ 
neous dressing G(/i,/l°,r) for 
AP transition at 802 nm 

^ Bphy 

25 nm 

half width half maximum ho¬ 
mogeneous dressing L(A, 4°, T) 
for BPh 


Table VI. Parameters for the simplified LHl model of R. rubrum. 


Parameter 

Value 

7 

10.06° 

Ay 

(Ay, -1- Ay 2 ) 

R 

(n -1- r2)/2 

d 

9D 

K 

2 

€ 

(fff + £/3)/2 

2 i> 22 

600, 377 cm-‘ 

„fJinhom 

36 cm“‘ 

std‘^^”' 

71 cm“‘ 

std”;"' 

277 cm-' 
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